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ABSTRACT 

Q \ We use the kinematics of satellite galaxies that orbit around the central galaxy in a 

■ dark matter halo to infer the scaling relations between halo mass and central galaxy 
' properties. Using galaxies from the Sloan Digital Sky Survey, we investigate the halo 
, mass— luminosity relation (MLR) and the halo mass— stellar mass relation (MSR) of 

central galaxies. In particular, we focus on the dependence of these scaling relations on 
the colour of the central galaxy. We find that red central galaxies on average occupy 
^ , more massive haloes than blue central galaxies of the same luminosity. However, at 

^■f-^ ' fixed stellar mass there is no appreciable difference in the average halo mass of red and 

\ blue centrals, especially for Af* ;< 10^°'^ /i^^M©. This indicates that stellar mass is a 

. better indicator of halo mass than luminosity. Nevertheless, we find that the scatter 

CO ' in halo masses at fixed stellar mass is non-negligible for both red and blue centrals. 

It increases as a function of stellar mass for red centrals but shows a fairly constant 

■ behaviour for blue centrals. We compare the scaling relations obtained in this paper 
' with results from other independent studies of satellite kinematics, with results from 

a SDSS galaxy group catalog, from galaxy-galaxy weak lensing measurements, and 
from subhalo abundance matching studies. Overall, these different techniques yield 
MLRs and MSRs in fairly good agreement with each other (typically within a factor 
of two) , indicating that we are converging on an accurate and reliable description of the 
\^ ' galaxy-dark matter connection. We briefly discuss some of the remaining discrepancies 

among the various methods. 

Key words: galaxies: halos — galaxies: kinematics and dynamics — galaxies: struc- 
ture — dark matter — methods: statistical 



1 INTRODUCTION 

The growth of structure in the Universe is predominantly 
driven by dark matter. The fluctuations in the dark mat- 
ter density fleld grow under the action of gravity and form 
a web-like structure. Galaxies form as baryon condensates 
at the density peaks of this cosmic web. Understanding the 
connection between the distribution of galaxies and the un- 
derlying distribution of dark matter is crucial to understand 
the physics of galaxy formation. This galaxy-dark matter 
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connection is often expressed in terms of the scaling rela- 
tions between the properties of galaxies and the mass of 
the dark matter halo in which they reside. Reliable mea- 
surements of the dark matter halo mass are essential to 
quantify these scaling relations. This can be accomplished 
with the help of numerous methods. These include tech- 
niques that are primaril y used on individu al systems such 
as rotation curves (e.g., iRubin et al. 1982t). stron g lensing 
of background galaxies (e.g., Gavazzi et al.ll2007|). and X- 
ray emision from hot gas in clust ers (e.g., Rvkoff et al.ll2008l : 
iDai. Kochanek. fc Morganll2007l ). With the advent of large 
scale galaxy redshift surveys, substantial progress has been 
made with methods that allow the inference of the dark mat- 
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ter halo masses in a statistical sense, e.g. the average halo 
mass as a function of various properties of galaxies. Such 
methods i nclude the modelling of the clustering of galax- 
ies ( e. g., Yang. Mo. fc van den Bosch 2003.: Zehavi et al.l 
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20051: iTinker et all l2005l: ICoUister fc Lahavl l2005l: 



Skibba et al]l2006l: Ivan den Bosch et al.ll2007l : iBrown et all 
20081: ISkibba fc She th| I2009I'). galaxy-galaxy weak lensing 



' '200d: 



Parker et al 



McKay et al. 2001 



Mandelbaum et al 



2007: .Mandelbaum. Seliak. fc Hirata 



20081 : ISchulz. Mandelbaum. fc PadmanabharJ 20091) and a 



combination of the t wo (e.g., Yoo et aLll2006l : Cacciato et al.l 
l2009l : lLi et al]|2009D . 

The satellite galaxies that orbit within the dark mat- 
ter haloes of their central galaxies are also excellent 
probes of the dark matter halo mass. Their kinemat- 
ics reflect the depth of the dark matter potential well 
they orbit. The number of satellites in massive systems 
like clusters is large enough to obtain a reliable mea- 
sure of their kinematics and hence the halo mas s (e.g., 
ICarlberg et al.l 19961 : ICMlberg. Yee fc Ellingsonlll997l '). How- 
ever, in low mass systems, where only a handful of satel- 
lites can be detected per central, one has to adopt a 
stacking procedure to quantify the kinematics of satellites 
llErickson. Gottesm an fc H unter 19871; Zaritskv et al.lll99j : 
IZaritskv fc Whitei l994: Zar itskv et'aLlll997l '). Central galax- 
ies with similar properties (e.g. luminosity) are stacked to- 
gether and the velocity information of their satellites is 
combined to obtain a quantitative measure of the kine- 
matics of the satellites. Various studies of the kinematics 
of satellite galaxies have now been carried out using large 
redshift surveys to improve the sample size of satellites 
llMcKav et al.ll2002| : iBrainerd fc Specianll2003l: iPrada et al ' 



20031 : iConroy et all l2005l: iBecker et al l l2007l : IConroy et al 
2OO7I : iNorberg. Frenk fc Cold 120081 ). The sample sizes ir 



these studies have been limited by their use of strict isola- 
tion criteria to identify centrals and satellites, specifically 
design ed to avoid misidentifications. Ivan den Bosch et al.l 
l|2004l ) devised relaxed selection criteria which were itera- 
tively adapted to the luminosity of central galax;ies to cir- 
cumvent this problem. Their criteria improved the sample 
size by nearly an order of magnitude over studies that use 
strict isolation criteria while still maintaining low levels of 
contamination. The improved statistics have warranted a 
better study of the systematics an d selection effects that 
bias the kinematic measurements (jNorberg. Frenk fc Cold 
I2OO8I ). 

In More, van den Bosch fc Cacciato (2009b; hereafter 
Paper I), we showed that if the relation between the halo 
mass and the stacking property has a non-negligible scatter 
then the kinematics of the satellites of the stacked system 
can be difficult to interpret. This issue has been neglected 
by most previous studies. We presented a new method to 
infer both the average halo mass and the scatter in halo 
masses as a function of the property used to stack the cen- 
tral galaxies. In More et al. (2009a; hereafter Paper II), 
this method was applied to galaxies from the Sloan Digital 
Sky Survey (York et al. 2000; hereafter SDSS) to infer the 
halo mass— luminosity relation of central galaxies (hereafter 
MLR). It was found that both the average and the scatter 
of the MLR of central galaxies increases with the luminosity 
of the central galaxy. 

The scatter in the MLR is an interesting quantity, as 



it is related to the stochasticity in galaxy formation. Two 
important, related questions are: (i) what is the physical 
origin of this stochasticity, and (ii) what galaxy property 
is most closely related to the mass of the halo in which it 
resides (i.e., shows the least amount of scatter at a given 
halo mass). The answer to (i) yields valuable insight into 
the physics of galaxy formation, while the answer to (ii) 
identifies the optimal galaxy property to trace the cosmic 
density field. It is well known that galaxies of the same stel- 
lar mass may have very different luminosities, even after 
correction for dust extinction. Galaxies with younger stellar 
populations will typically be bluer and more luminous than 
galaxies of the same stellar mass, but with an older stellar 
population. It may well be that the stellar mass of a central 
galaxy is a better indicator of halo mass than its luminosity, 
in which case the halo mass— stellar mass relation of cen- 
tral galaxies (hereafter MSR) will have less scatter than the 
MLR, and the scatter in the MLR will be correlated with 
the color of the central galaxy. Obviously, since we lack a 
complete theory of galaxy formation, it may also be that 
the opposite holds, and that the scatter in the MLR is actu- 
ally less than that in the MSR. In this paper, we investigate 
these issue by measuring the kinematics of satellite galaxies 
as functions of both the luminosity and stellar mass of cen- 
trals split by colour into red and blue sub-samples. Using 
the methodology outlined in Paper I, we use these to probe 
both the means and scatters of the MLR and MSR of red 
and blue galaxies. 

This paper is organized as follows. In Section (2] we de- 
scribe the data used in this paper. In Section |3] we explain 
our method of analysis. In particular, we describe the pro- 
cedure used to identify the centrals and the satellites, the 
measurement of the kinematics of satellites, and the sub- 
sequent modelling to determine the halo masses of central 
galaxies. In Section |4] we present our results and compare 
them with other independent studies. Finally, we summa- 
rize our findings in Section (5) 

Throughout this paper we adopt the cosmological pa- 
rame ters supported by the 3 year data release of WMAP 
(Spe rgel et all l2007h : Q.^ = 0.238, = 0.762, h = 
i?o/100 kms~^ Mpc~^ = 0.734, the spectral index of ini- 
tial density fluctuations Us = 0.951 and the normalization 
of the power spectrum of density fluctuations erg = 0.744. 
We use the symbol M to refer to the mass of a dark matter 
halo, which is defined as the mass enclosed within a spherical 
overdensity 5p/p = 200, where p denotes the mean matter 
density of the universe. 



2 DATA 

We use data from the SDSS which is a joint flve-passband 
{u, g, r, i and z) imagi ng and medium re solution (R ~ 1800) 
spectroscopic survey (|York et al.l I2OO0I ). More speciflcally, 
we use the New York Univers ity Value Added Galaxy Cat- 
alogue (|Blanton et al.l 2005bll. which is based up on SDSS 
Data Release 4 l|Adelman-McCarthy et al.l 120061 ) but in- 
cludes a set of significant improvements over the original 
pipelines. The magnitudes and colours of the galaxies are 
based upon the standard SDSS Petrosian technique and have 
been k-corrected and ev olution corrected to z = 0.1 using 
the method described in iBlanton et ai] (|2003al lb[). The no- 
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tations '^'^(g — r) and "'^Af,. — 51og/i are used to denote 
the resulting {g — r) colour and the absolute magnitude of 
the galaxies. From this catalogue, we select all galaxies in 
the main galaxy sample with apparent magnitudes less than 
17.77 that lie in an area where the redshift completeness 
limit of the survey C > 0.8. Next we construct a volume lim- 
ited sample that is complete in luminosity above a "'^r-band 
luminosity of Lmin = 10^'^ /i~^L(rl3. The redshift range that 
we adopt for this volume limited sample is 0.02 ^ z ^ 0.072, 
which results in a total sample of 58, 396 galaxies. 

Stellar masses are indicated by M« and are computed 
using the relation between the ste llar mass-t o -light ratio and 
the ° °{g - r) colour provided bv lBell et all l|2003l ): 



log 



fe-2M0 



-0.306 + 1.097 [°-'\g - r)] - 0.10 
-0.4(''°A/r-51og/i-4.64). (1) 



Here ""{g — r) and '''"Mr — 5 log h denote the (g — r) colour 
and the r-band absolute magnitude of galaxies k-corrected 
and evolution corrected to z = 0.0, 4.64 is the r-band mag- 
nitude of the Sun in th e AB system, a nd the —0.10 term is 
a re sult of ad opting the Kro upal (|200lh initial mass function 
(see lBorch et a l. 2006). The typical uncertainties in the stel- 
lar masses obtained i n this manner are of the order of ~ 0.1 
dex (jBell et al.ll200^ . 

We classify galaxies to be red or blue based upon their 
bimodal distribution in the ""(g — r) colour-stellar mass 
plane. We use the following separation criterion to demar- 
cate the boundary between red and blue galaxies in the 
colour-stellar mass plane (see Appendix E)) : 

"■"(5 - r)cut = 0.65 + 0.10 (log [M,/{h-^MQ)] - lO.O) .(2) 



3 METHODOLOGY 
3.1 Selection criteria 

The first step towards measuring the kinematics of satellite 
galaxies in the SDSS is to accurately identify central galaxies 
and their associated satellites. For our analysis of the MLR 
we proceed as follows. A galaxy is identified to be a central 
if it is brighter than every other galaxy within a cylindrical 
volume specified by ii < Rh and \AV\ < (AV)h centred on 
itself. Here R is the physical distance from the galaxy under 
consideration projected on the sky and AV is the line-of- 
sight (hereafter los) velocity difference between two galaxies. 
All galaxies that lie within a cylindrical volume specified 
by i? < -Rs and lAV'l < {AV)s around a central galaxy, 
and that are fainter than the central galaxy, are labelled to 
be its satellites. The criteria used to select the sample of 
central and satellite galaxies for the analysis of the MSR are 
almost identical, except that in this case the central galaxy 
must have the largest stellar mass in its cylindrical volume 
specified by Rh and {AV)h- 

The parameters Rh,{AV)h,Rs and {AV)s define 
the sizes of the cylinders used to identify central 



galaxies and their satellites. Contr ary to most previ- 
ous studies of satellite kinematics (iMcKay et al.l I2OO2I: 
'Brainerd & Specian 2003; 'Prada ct al. 20031: Iconrov et al l 
.2005, , 2007 ; Norberg , Frcnk fc Cole 2008), we do not use 
fixed values for these parameters. Rather, since halo mass 
is expected to be positively correlated with the luminos- 
ity or stellar mass of the central galaxies, we scale the se- 
lection parameters according t o the property of the galaxy 
under consideration. Following Ivan den Bosch et al.1 (|2004l ). 
we adopt Rh = 0.80-200/1"^ Mpc, {AV)h = 1000cr2oo kms"\ 
Rs = Q.15a2ooh-^ Mpc and (AF)s = 4000 kms'^ Here 0200 
is the satellite velocity dispersion in units of 200 kms~^, 
which we parameterize as 



o-20o(log Qio) = a + b (log Qio) + c (log Qio) 



(3) 



where Qio is either the central galaxy luminosity in units 
of 10^° h~^LQ or the stellar mass in units of lO^°/i~^M0, 
depending upon the property used to stack central galaxies. 
Clearly, since the determination of cr2oo requires a sample 
of centrals and satellites, this selection method has to be 
iterative. Fixed values of the selection criteria parameters 
are used to identify the central and the satellite galaxies in 
the first iteration. The velocity dispersion of the selected 
satellites as a function of the central galaxy property, pa- 
rameterized via Eq. (jSJ, is fit using a maximum likelihood 
method and subsequently used to scale the values of the pa- 
rameters that define the selection criteria. These are used to 
select a new sample of centrals and satellites, and the en- 
tire procedure is repeated until convcrgenc fl Using detailed 
mock galaxy redshift surveys, [van don Bo sch et al.l (|2004l ) 
have shown that this iterative technique yields much lower 
interloper fractions than the more common method using 
fixed cylindrical volumes (see also Paper II). For complete- 
ness, and to allow the reader to reproduce our results. Ta- 
ble 1 lists the final iteration criteria used for our various 
samples (in terms of the parameters a, b, and c that ap- 
pear in Eq. [3]), as well as the total number of centrals and 
satellites selected in each sample. Note that these parame- 
ters differ depending on the centrals we choose to stack for 
our analysis. The sample of centrals and satellites selected 
when the selection criteria are tuned based on the velocity 
dispersion around all centrals stacked by luminosity (stellar 
mass) is called Sample LA (SA). Samples LR (SR) and LB 
(SB) are selected by tuning the selection criteria parame- 
ters based upon the velocity dispersion around red and blue 
centrals stacked by luminosity (stellar mass), respectively. 



3.2 Velocity dispersion measurement 

In Paper I we demonstrated that the commonly measured 
velocity dispersion of satellite galaxies, Osat, cannot be used 
to uniquely determine the scaling relation between halo mass 
and a central galaxy property unless this relation has zero 
scatter. In fact, we have shown that different scaling rela- 
tions with different amounts of scatter can yield exactly the 
same Osax- In the same paper, however, we have shown that 
this degeneracy can be broken using a combination of two 
different measures for the velocity dispersion of the satellite 



^ Th e '-'•^ r-band magnitud e of the Sun in the AB system equal to 
4.76 l|BIanton et al 1 l2003sJ ) is used to convert the absolute mag- 
nitude of a galaxy to its luminosity in units of h~^'LQ. 



We refer the reader to Paper II and Ivan den Bosch et al.l ||2004| ') 
for details regarding this method. 
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Table 1. Selection criteria parameters 



Samples 


a 


b 


c 


Centrals 


Satellites 


LA 


2.19 


0.37 


0.30 


3949 


6213 


LR 


2.23 


0.34 


0.31 


2723 


4873 


LB 


2.11 


0.47 


-0.08 


1082 


1255 


SA 


2.07 


0.22 


0.20 


3834 


6232 


SR 


2.11 


0.19 


0.20 


3095 


5445 


SB 


1.99 


0.49 


-0.20 


701 


786 



The parameters a, b and c that define the criteria used to select 
central and satellite galaxies for all the samples used in this paper 
(see text), and the total numbers of centrals and satellites thus 
selected. 



galaxies: satellite- weighted (cTsw) and host- weighted (cr^^). 
To measure the velocity dispersion in these two different 
schemes, the satellite galaxies in the final sample are first 
binned into sub-samples based upon the properties (lumi- 
nosity or stellar mass) of their central galaxies. For each 
bin, the distribution of los velocities of satellite galaxies with 
respect to their centrals, P{/S.V), is constructed by either 
giving each satellite equal weight (satellite-weighting) or a 
weight equal to 1/A^sat (host-weighting), where A^sat denotes 
the number of satellites around the host of the satellite under 
consideration. As shown in Paper I, the difference between 
iTsTO and cr^w depends on the amount of scatter in the scal- 
ing relation between halo mass and central galaxy propertjjf], 
and allows the degeneracy to be broken. 

In order to extract the satellite velocity dispersion 
(satellite-weighted or host-weighted) from the correspond- 
ing P(A1/) distributions, we fit P{AV) using the sum of 
two Gaussians plus a constant: 



P(A1/) = ao + ai exp 



2(7? 



+ a2 exp 



2ol 



•(4) 



The satellite velocity dispersion (satellite-weighted or host- 
weighted) then follows from 



(sw/hw) 



a\ CTi -I- a2 (72 



(5) 



Here tJcrr is the contribution to the effective variance of 
P(A1/) due to redshift errors in the SDSS. Given that each 
individual galaxy has a redshift error of ~ 35kms~^, the er- 
ror on the velocity difference, AV^, of the central and satel- 
lite galaxies IS fJ err — a/2 X 35 ~ 49.5 kms \ which is 
the value we adopt throughout. The errorbars on the veloc- 
ity dispersions were estimated as the variance in the veloc- 
ity dispersions obtained by fitting 1000 realisations of the 
P(AF) histogram (generated assuming Poisson errorbars) 
with the procedure described above. Detailed tests using 
mock catalogues have shown that the above method yields 
extremely reliable estimates of the actual velocity disper- 
sions (see Paper II). 

Applying this method to the central-satellite samples 
selected using the criteria described in Section 13.11 we ob- 
tained the satellite-weighted and host-weighted velocity dis- 
persions as well as the average number of satellites per cen- 
tral as a function of luminosity for samples LA, LR and LB, 



and as a function of stellar mass for samples SA, SR and 
SB. In addition, we also measured the fraction of red cen- 
tral galaxies as a function of luminosity from Sample LA 
and as a function of stellar mass from Sample SA. In the 
following subsection we describe how to use these data to 
constrain the MLR and MSR of central galax;ies. 

3.3 The Model 

The satellite-weighted and host-weighted velocity disper- 
sions of satellites and the average number of satellites of 
a central galaxy with a given property depend on the distri- 
bution of halo masses of central galaxies with that property, 
P(M\Q'). Therefore these observables can be used to infer 
the mean and the scatter of the scaling relation between 
halo mass and the central galaxy property. The analytical 
expressions that describe these three quantities are given by 
(see Paper I): 



XQ) = 



\^ P{M\Q) {iVe)ap,M (aDap.M dM 



2 

0"hv 



(Q) = 



/■ex 




P{M\Q) (Afs)ap,AfdM 
P{M\Q)V{M) (aDap.MdM 



{Ns){Q) 



P{M\Q)r{M)dM 
/^°°P(M|Q)(Ar,)ap.AfdM 



PiM\Q)r{M)dM 



(6) 



(7) 



(8) 



Here {A's)ap,A/ and {(Tsat)ap,A/ denote the average number 
of satellites and the average velocity dispersion of satellites 
within the aperture in a halo of mass M, respectively. The 
factor 'P(M) is the fraction of haloes of mass M that host at 
least one satellite. Our way of modelling the observables thus 
consists of two parts: (i) the kinematics of satellite galaxies 
in a single halo of a given mass and (ii) the statistics that 
describe how central and satellite galaxies occupy haloes. 
We describe each of these parts in the following subsections. 



3.3.1 Kinematics in a single halo 

We assume that dark matter haloes are spherically symmet- 
ric and that t heir density distribut ions follow the universal 
NFW profile ()Navarro. Frenk. fc White„1997. ). 



p{r\M) oc 



(9) 



where Vs is a scale radius specified in terms of the virial 
radius rvir using the concentration parameter , c = ryir/r^. 
We us e the concentration-mass relation from iMaccio et al.l 
(j2007h . appropriately modified for our definition of the halo 
mass. We assume that the number density distribution of 
satellite galaxies, nsat(r|M), is given by 



".sat(riM) oc 



f— 



1 + 



r 



7-3 



(10) 



^ In the case of zero scatter, one has that tTs, 



Here 7 represents the slope of the number density distri- 
bution of satellites as r — > and 7?. is a free parameter. 
Throughout this paper, we use the result from Paper II that 
the number density distribution of satellites in the SDSS 
can be well described by Eg. ifTOll with 7 = 0.0 and TZ = 2 
(see also lYang eral]|2005l and lMorell20093 ). The distribution 
nsiit{r\M) is normalized such that 
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{Ns)(M)=47T 



nsat(r|M) dr , 



(11) 



where {Ns){M) denotes the average number of sateUites in 
a halo of mass M. The number of sateUites within the aper- 
ture, {Ns)!!,p,M, is then given by 



{Ns) 



a.p,M 



= 4tv 



RdR 



nsat(r|M) 



rdr 



(12) 



/O JR V^^^W 

We further assume that the satellite occupation numbers 
follow Poisson st atistics, which is supporte d both by direct 
observations (e.g. lYang. Mo fc van den Bo sch 2008) and by 
numerical simulations i Kravtsov et al.l 2004 ). The fraction 
of central galaxies that have at least one satellite within the 
aperture radius is then given by 



V{M) 



- exp[-{7Vs)ap,A/] 



(13) 



To find an analytical expression for {cTss.t)a.p,M, first note 
that the Jeans equation can be used to find an expression for 
the radial velocity dispersion of satellites at a radial distance 
r from the centre; 



f — 



1 + 



3-7 



Here Vvir is the circular velocity at r^-^ and 



(14) 



/ y(l + y)~'dj/, (15) 

"'0 

(see Paper I for a detailed derivation). If the velocity dis- 
persion of the satellites is assumed to be isotropic, then 
((Tgat)ap,M cau be expressed as the average of the radial ve- 
locity dispersion, o'^^^(r\M)^ over the aperture R^: 

RAR 



{''"sat)ap,M — 



47r 



rdr 



I ™ n.at(r|Af) aL(r-IM) . 



(16) 



Note that the denominator has no dependence on M and is 
just a multiplicative normalisation constant which cancels 
out when we model our observables (see Eqs. [S]-[S1). 

We use simple parametric forms to model the distribu- 
tion P{L\M^C) and /c(M). The distribution P(L|AfnC) 
is assumed to be a log-normal given by 



V27rcr log L 



- exp 



(iog[^/£])^ 



dL 
T 



(18) 



Here log L(C,M) denotes the mean of the log-normal dis- 
tribution and aiog l (C) is the corresponding scatter. We use 
four parameters each to specify the relations L{R, M) and 
L{B,M): a low mass end slope, 71, a high mass end slope, 
72, a characteristic mass scale, Mi, and a normalisation, Lq; 



{M/Mi 



[l-f (M/Mi)]^i-^2 ■ 



(19) 



We further assume that the scatters aiog l{R) and aiogL(B) 
are independent of halo mass. Thus for each colour C, we 
use 5 parameters to describe the distribution P{L\M n 
C). This parametrization i s mot ivated by the results of 
lYang, Mo fc van den BoschI (|2008l '). who inferred the con- 
ditional luminosit y function using t he large SDSS galaxy 
group catalogue of lYang et al.l (|2007| ). 

The function /_r(M) is assumed to be linear in logM: 

f'n{M) = fo + af[log{M/h-'MQ)-12.0] 

faiM) = min(max[0,/;i(M)],1.0) (20) 

where the second equality takes into account that fn is 
a fraction, and therefore bounded by zero and unity (cf. 
Ivan den Bosch et all [20031) . Also note that /b(M) = 1 - 
/fl(M), as 'red' and 'blue' form a mutually exclusive and 
exhaustive set of colours assigned to the central galaxies. 
Hence, /_r(M) and /b(M) add a total of two free parame- 
ters to our model. 

We also need a model for the satellite occupation num- 
bers, {Ns}{M). Throughout we assume that the number of 
satellite galaxies in a halo of mass M scales with halo mass 



3.3.2 Halo occupation statistics 

To model the kinematics of satellite galaxies around central 
galaxies stacked by a particular property, we need to specify 
the distribution of halo masses for central galaxies with that 
property, P{M\Q). Let us first consider the case of central 
galaxies stacked by luminosity and colour. In this case, Q = 
L n C, where L denotes the luminosity of the central galaxy 
and C its colour. We will use the letters R for red and B 
for blue when referring to the actual colours. Using Bayes' 
theorem. 



P{M\L n C) 



P{L\MnC)fc{M)P{M) 



(17) 



fciL)P{L) 

Here P{L\Mr]C) describes the distribution of central galaxy 
luminosities in haloes of mass M that host central galaxies 
of a particular colour C, fc{M) = P{C\M) is the fraction 
of haloes of mass M that host a central galaxy of colour C, 
and P(M) is proportional to the halo mass function n(Mj3. 



(iVs)(A/) = iVi2 



M 



1012 h-^Mr. 



(21) 



which adds two more parameters to our model; N12 and 
a. Note that we assume that {Ns}{M) is independent of the 
colour of the central galaxies. Although we believe this to be 
a realistic assumption, we will discuss the potential impact 
of its violations in Section [l] 

Next, consider the case where galaxies are stacked only 
according to their luminosity, i.e., Q = L. To model the 
observables we need to know the distribution 



PiMlL)^mM)IiMl^P^^M)niM). 



(22) 



The distribution P{L\M) is related to the distributions 
P(L n R\M) and P{L n B\M) according to: 

P(L|A'/) = P(L|Af n R)fR{M) + P(LjAf n B)fB{M) (23) 



^ For the analysis i n this paper, we use the halo mass function of 
[Tinker et al. I liooi) for which haloes are defined as spheres with 



an average density that is 200 times the average matter density 
in the universe. 
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Finally, the expression that describes the fraction of red cen- 
trals as a function of central galaxy luminosity is given by 



fB{M) S{A'L,B) P{AU\M n B) 



(27) 



P{L\M n 7?) V{M) n{M) fR{M)dM 
P(L\M) P{M) n{M)dM 



(24) 



Note that we have appropriately corrected for the fact that 
the observed fraction of red centrals is calculated using only 
those centrals that have at least one satellite. 

Hence, our analytical model has a total of 14 free pa- 
rameters and completely describes the kinematics of satellite 
galaxies and the average number of satellites around centrals 
stacked by luminosity in Samples LA, LR and LB. This an- 
alytical framework also allows us to calculate the fraction of 
red centrals as a function of luminosity. 

For the analysis of the kinematics of satellite galaxies 
around centrals stacked by stellar mass, there is an addi- 
tional complication that has to be addressed. The central 
and satellite galaxies used for our analysis are selected from 
a volume limited sample that is complete above a certain 
luminosity. Since we have used both colour and luminosity 
to assign the stellar masses, our sample starts to become 
incomplete in stellar mass roug hly below 10^°/i"^Mq. The 
completeness is a function of both stellar mass and colour, 
and is described by the sample selection function S{Mi,, C), 
defined as the fraction of galaxies of stellar mass Af, and 
colour C in the SDSS volume with 0.02 2 ^ 0.072 (i.e., 
our sample volume) that make it into the sample. The de- 
termination of S{M,,C) is discussed in Appendix El 

The sample selection function ^(Mt , C) enters our 
model in the following way. We can write the stacking prop- 
erty as Q = M, nCnS, where we use S to denote the subset 
of all galaxies with stellar mass A/* and colour C in our sam- 
ple volume that make it into the sample. The corresponding 
distribution of halo masses, P(Af jAf, nC C] S), can then be 
written as 



P(A^|Af, n cn S) = 

_ g(Af.,C) 



P{S\M n Af, n C) P{M n Af, n C) 



P(Af. n c n S) 



P{M, n c n 5) 

fciM) P(M) P{M,\M n C) . 



(25) 



In the second equality, we have identified the distribution 
P{S\M n Af, n C) as the selection function, S{M^,C). In 
addition, we have also expressed P{M D n C) in terms 
of P{M,\M n C). Note that the selection function does not 
depend on the halo mass M and acts as a multiplicative nor- 
malisation constant for the distribution P{M\M,C\CC\S). It 
harmlessly cancels out from the expressions that analytically 
describe the observables when central galaxies are stacked 
by colour. However, it turns out to be important for calcu- 
lating P{M\Mi, nS) and hence the observables when central 
galaxies are stacked by stellar mass alone. First note that 
the distribution P{M\M, n S) is related to P(A'/.|A^ n S) 
such that 



P{M\M, r\S) = 



P{M, nS\M)P{M) 
P{M, n S) 



(26) 



The distribution P{Mt n S\M) can be expressed in terms of 
P{Mt\M n R) and P{M,\M n B) as follows: 

P{M, n S\M) = P{M, nRn S\M) + P{Ah nBn S\M) 

= /fl(Af) S{M,,R) P{KU\M nR) + 



Similar to Eq. (fT5)) . we parameterize P(Af, |A/ n C) as 
a log-normal distribution with mean log Af, (C, A4") and a 
scatter aiogAf, which depends on colour but is indepen- 
dent of halo mass. The relation A/«(C, M) is described using 
Eq. (|19|l , but with L replaced by Mt , and we assume that the 
scatters criogA/, (P) and criogA/. (B) are independent of halo 
mass. For fii{M) and {Ns){M) we adopt the same param- 
eterizations as before (i.e., Eqs. [50] and [H]). Hence, our 
model for the analysis of the MSR also contains 14 free pa- 
rameters, which we constrain using the observables obtained 
by stacking central galaxies by stellar mass and colour. 



3.4 Constraining the model parameters 

We now describe our method to constrain the model param- 
eters. Here we focus on the analysis of the MLR, but note 
that the analysis of the MSR is basically the same. We have 
measurements of the satellite- weighted velocity dispersion, 
the host-weighted velocity dispersion and the average num- 
ber of satellites per central, for 10 different luminosity bins, 
and for each of the three samples LA, LR and LB. In ad- 
dition, we have 10 measurements of the fraction of red cen- 
trals as a function of luminosity. Since most of the centrals 
in sample LA are present in either sample LR or sample 
LB, the velocity dispersions and average number of satel- 
lites measured from sample LA are not independent from 
those obtained using samples LR and LB. Therefore, we do 
not use these measurements from sample LA to constrain 
the model parameters. This leaves a total of 70 independent 
data points to constrain our 14 model parameters. 

We use flat uninformative priors on each of the model 
parameters (albeit in a limited interval for each of the pa- 
rameters). We use a Monte-Carlo Markov Chain (hereafter 
MCMC) technique to sample from the posterior probability 
distribution of each of these parameters given the observa- 
tional constraints. The MCMC is a chain of models, each 
consisting of the 14 parameters. At any point in the chain, a 
trial model is generated with the 14 free parameters drawn 
from a 14-dimensional Gaussian proposal distribution which 
is centered on the current values of the parameters. The 
chi-squared statistic, Xtry, for this trial model, is calculated 
using 



2 

Xtry 



J2 [xUC) + Xhw(C) + xUc)] + xl 



C=R,B 



where 

10 

Xsw(C) = J2 
i=l 

10 r 

Xhw(C) = Yl 

1=1 

10 

xUc) = 

10 

E 



a-sw(Q[i]) — d-sw{Q[i]) 



2 

Xb 



A<Tsw((3[i]) 

o'hw(QH) - o-hw(Q[i]) 
A<Thw(QM) 

{iVs)(QH)-iVs(QH) ' 
ANs{Q\i]) 



fR{L)-fR{L\i]) 



A/b(LW) 



(28) 

(29) 
(30) 
(31) 
(32) 



Here Q = L C] C , X denotes the observable X and AX 



Table 2. MLR: Percentiles of the posterior distributions 
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Table 3. MSR: Percentiles of the posterior distribution 



Parameter 16 percent 50 percent 84 percent 



Red log(Lo) 9.62 9.87 10.24 

centrals log(A/i) 11.20 11.55 11.96 

71 2.42 3.37 4.43 

72 0.26 0.36 0.43 
(TiogL 0.18 0.20 0.23 



Blue log(Lo) 8.98 9.32 9.67 

centrals log(Afi) 10.14 10.45 10.91 

71 2.29 3.08 4.26 

72 0.32 0.46 0.58 
o-iogL 0.13 0.20 0.26 



Parameter 16 percent 50 percent 84 percent 



Red log(Af.o) 10.12 10.63 10.97 

centrals log(Mi) 11.42 11.94 12.29 

71 2.44 3.44 4.49 

72 0.15 0.29 0.44 
o-logM. 0.16 0.19 0.22 



Blue log(M*o) 8.36 9.42 10.68 

centrals log(Mi) 10.57 11.29 11.98 

71 2.32 3.21 4.37 

72 0.48 0.98 1.31 
o-iogM. 0.08 0.15 0.27 



0.37 0.45 0.56 

0.43 0.61 0.77 

-0.91 -0.79 -0.68 

1.18 1.27 1.38 



The 16, 50 and 84 percentile values of the posterior distributions 
for the parameters of our model obtained from the MCMC anal- 
ysis of the velocity dispersion data from Samples SR, SB and 
SA. 



for the average number of satellites. The fraction of red cen- 
trals as a function of luminosity obtained from Sample LA 
are shown in panel (g). The data in these 7 panels is used 
to constrain the 14 parameters of our model that describe 
the halo occupation statistics of red and blue centrals, and 
the satellite occupation numbers in haloes that host red and 
blue centrals. The blue and purple shaded regions indicate 
the 68 and 95 percent confidence levels obtained from our 
MCMC. A comparison with the data reveals that the model 
accurately reproduces the data. 

Fig. [2] shows the velocity dispersions and the average 
number of satellites as a function of luminosity around all 
centrals obtained from Sample LA. As for samples LR and 
LB, the model is in excellent agreement with the data (open 
squares), even though these data were not directly used to 
constrain the model. 

The 16, 50 and 84 percentile values of the posterior dis- 
tribution for our model parameters obtained from our analy- 
sis are listed in Tabled The constraints on our model ingre- 
dients are presented pictorially in Fig. [31 Panels (a) and (d) 
show the constraints obtained on log L{M) for red and blue 
centrals, respectively, while panels (b) and (e) show the cor- 
responding posterior distributions of the scatter in luminosi- 
ties at fixed halo masses. Panel (c) shows the average number 
of satellite galaxies as a function of halo mass, and panel (f) 
shows the fraction of red centrals as a function of halo mass. 
At the bright end, the mean luminosity of red central galax- 

36+° '"' 

ies scales with halo mass as L (x M ' while that of 

blue central galaxies scales as L cx 

M^-^^-o.n, At the faint 
end, the constraints on the faint end slope of the logI/(Af) 
relation are entirely dominated by the prior 71 G [2.0,5.0]. 
This is due to the magnitude limit of the SDSS, which 
is not sufficiently faint to reliably probe the occupation 
stat istics of dark matter haloes with M <, 10^^ M© (but 
see lYang. Mo fc van den BoschI |2009| ). For the scatter in 
P{L\M) we obtain that criogL = 0.20lg;°2 ^or red centrals 
and (Tiog L ~ 0.20lIiQ g7 for blue centrals. Note that the scatter 



All fo 0.57 0.64 0.71 All /o 

centrals af 0.14 0.20 0.27 centrals 

logAfi2 -0.83 -0.72 -0.62 logAfi2 

a 1.19 1.28 1.38 a 



The 16, 50 and 84 percentile values of the posterior distributions 
for the parameters of our model obtained from the MCMC anal- 
ysis of the velocity dispersion data from Samples LR, LB and 
LA. 

its corresponding error. The trial step is accepted with a 
probability given by 

p J 1.0^ if Xtry ^ Xcur /'QQ'l 

/^accept - I exp[-(x?.y - X?u.)/2], if Xly > xLr ^ ' 

where Xcur denotes the x^ for the current model in the chain. 

We initialize the chain from a random position in our 
14-dimensional parameter space and discard the first 20, 000 
models (the 'burn-in' period) allowing the chain to sample 
from a more probable part of the distribution. We proceed 
and construct a chain consisting of 40 million models. We 
thin this chain by a factor of 10"^ to remove the correla- 
tions between neighbouring models. This leaves us with a 
chain of 40, 000 independent models that sample the poste- 
rior distribution. We use this chain of models to estimate 
the confidence levels on the parameters and the relations 
of interest, namely the mean and the scatter of the scaling 
relation between halo mass and the central galaxy property 
under consideration. 



4 RESULTS 

4.1 The Halo Mass— Luminosity Relation 

The analysis of the MLR of central galaxies is carried out 
by analyzing Samples LA, LR and LB (see Table 1 for the 
selection criteria and the numbers of centrals and satellites 
in each of these samples). The host-weighted velocity dis- 
persion, the satellite-weighted velocity dispersion and the 
average number of satellites as a function of the luminos- 
ity of the central galaxies obtained from Samples LR are 
shown as open squares in panels (a), (b) and (c) of Fig. \T\ 
respectively. In the same figure, panels (d), (e) and (f) show 
the corresponding observables obtained from Sample LB. At 
fixed luminosity, the velocity dispersion of satellite galaxies 
around red centrals are systematically larger than the veloc- 
ity dispersions around blue centrals. The same is also true 
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Figure 1. Obscrvablcs used to constrain the MLR of central galaxies (open squares with crrorbars). The upper and middle panels show 
the observables measured using red centrals (Sample LR) and blue centrals (Sample LB), respectively. From the left to the right these 
panels show the host- weighted velocity dispersions [panels (a) and (d)], the satellite- weighted velocity dispersions [panels (b) and (e)], 
and the average number of satellites per central [panels (c) and (f)], all as function of the luminosity of the central. Panel (g) shows the 
fraction of red centrals as a function of luminosity as measured from Sample LA. The blue and purple regions indicate the 68 and 95 
percent confidence intervals obtained from the MCMC, showing that the model accurately fits the data. 




log (Lyh-^L,,) log (Lyh-^L^) log (Lyh-^L,,) 

Figure 2. Obervables measured using Sample LA. Panels (a), (b) and (c) shows the host-weighted velocity dispersions, the satellite- 
weighted velocity dispersions, and the average number of satellites per central, respectively, all as a function of the luminosity of the 
central galaxy. Although these observables are not used to constrain the model parameters (see text), the 68 and 95 percent confidence 
intervals obtained from the MCMC, indicated by the blue and purple regions, show that the model accurately fits these data as well. 
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Figure 3. Constraints on tiie model ingredients obtained from thie MCMC. The 68 and 95 pereent confidence intervals are indicated as 
blue and purple regions, respectively. Panel (a) shows the average luminosity of red centrals as a function of halo mass, wrhile panel (b) 
shows the posterior distribution of the scatter cri^g l in this relation. Panels (d) and (e) show the same but for blue centrals. Panel (c) 
shows the constraints on the average number of satellites as a function of halo mass, and finally, panel (f ) shows the constraints on the 
fraction of red centrals as a function of halo mass. 
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Figure 4. Constraints on the MLR obtained from the MCMC. The average halo masses as a function of the luminosity of red centrals, 
blue centrals, and all centrals are shown in panels (a), (b) and (c) respectively. Panels (d), (e) and (f) show the confidence intervals on 
the scatter in these relations. As in the previous figures, blue and purple colours indicate the 68 and 95 percent confidence intervals. 
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Figure 5. Comparison between the average MLR (left hand panel) and the average MSR (right hand panel) of all centrals (squares), 
red centrals (triangles) and blue centrals (circles). Errorbars indicate the 68 percent confidence intervals. See text for discussion. 



for blue centrals is less well contrained than for red centrals, 
which is due to the smaller sample size (see Table 1). 

Fig. |4] shows the average halo mass and the scatter in 
halo masses obtained from our analysis for red centrals (left- 
hand panels), blue centrals (middle panels) and all centrals 
(right-hand panels) as a function of luminosity. As expected, 
brighter galaxies reside on average in more massive haloes. 
Note that the scatter in halo masses around all centrals ap- 
pears somewhat higher than that around either red or blue 
centrals. This indicates that some fraction of this scatter is 
due to the fact that red and blue centrals of the same lumi- 
nosity reside, on average, in haloes of different mass (see be- 
low). However, red and blue centrals separately still reveal a 
significant amount of scatter in their halo masses, whereby 
the scatter around red centrals increases with luminosity, 
while that around blue centrals shows no clear luminosity 
dependence. 

The left-hand panel of Fig. [5] summarizes these results. 
It compares the MLR of red centrals (triangles) to that of 
blue centrals (circles) and to that of all centrals (squares). 
At fixed luminosity, red centrals on average reside in more 
massive haloes than blue centrals. As expected, the MLR of 
all centrals lies in between that of red and blue centrals. As 
the fraction of red centrals steadily increases to unity at the 
bright end, the average halo mass of all centrals shifts from 
tracing the MLR of blue centrals to tracing the MLR of red 
centrals. 



4.2 The Halo Mass Stellar Mass Relation 

The analysis of the MSR of central galaxies is carried out 
by analysing Samples SA, SR and SB (see Table 1 for the 
selection criteria and the numbers of centrals and satellites 
in each of these samples) . The host-weighted velocity disper- 
sion, the satellite-weighted velocity dispersion and the aver- 
age number of satellites as a function of the stellar mass of 
the central galaxies obtained from Samples SR are shown 
as open squares in panels (a), (b) and (c) of Fig. [S] re- 
spectively. In the same figure, panels (d), (e) and (f) show 
the corresponding observables obtained from Sample SB. At 



fixed stellar mass, the velocity dispersion of satellite galaxies 
around red centrals is systematically larger than that around 
blue centrals. The same is also true for the average number of 
satellites. The fraction of red centrals as a function of stellar 
mass obtained from Sample SA are shown in panel (g). The 
sharp drop of /r to zero at the low stellar mass end is due to 
our use of a volume limited sample complete in luminosity, 
which causes the sample selection function for red centrals 
to go to zero at the low mass end (see Appendix [XJ . The 
data in these 7 panels is used to constrain the 14 parameters 
of our model. The blue and purple shaded regions indicate 
the 68 and 95 percent confidence levels obtained from our 
MCMC. As in the case of the MLR, the model accurately 
fits the data, also in the case of sample SA, which is not 
used to constrain the model (see Fig. (Tjl. The 16, 50 and 84 
percentile values of our model parameters obtained from our 
analysis are listed in Tabled 

Fig. [8] is the same as Fig. |3] except that it shows the 
constraints on the MSR, rather than the MLR. At the mas- 
sive end, the average stellar mass of red centrals scales with 
halo mass as M, oc M ' -014 while for blue centrals we 

a+o.3.3 

For the scatter in P(A'/«|M) 



find that M. oc M°-®*-o ™. 
we obtain that oia^M, = 0.19~ 



C" log Mt 



= 0.15 



+0.12 



0.03 



for red centrals and 
for blue centrals. The posterior distri- 
bution of CTiogAf, for blue centrals has a long tail extend- 
ing to large values of scatter. At the low end, the confi- 
dence levels for aiogM, are largely dominated by our prior 
(ciogA/, > 0.04), which has been adopted for computational 
convenience. Hence, we cannot rule out that the data is con- 
sistent with zero scatter (i.e. (TiogAf* = 0.0) for blue centrals 
stacked by stellar mass. 

In Fig. [^we present the average halo mass and the scat- 
ter in halo masses obtained from our analysis for red centrals 
(left-hand panels), blue centrals (middle panels) and all cen- 
trals (right-hand panels) as a function of stellar mass. Note 
that the scatter in halo masses increases with stellar mass 
for the red centrals, while it is roughly independent of stellar 
mass for blue centrals. In the right-hand panel of Fig. [S] we 
compare the MSR of red centrals (triangles) to that of blue 
centrals (circles) and to that of all centrals (squares). At 
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Figure 6. Same as Fig. [T] but as function of the stellar mass of the centrals, based on samples SR, SB and SA. 
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Figure 7. Same as Fig. [2] but as function of the stellar mass of the centrals, based on sample SA. As for Fig. [2] these observables are 
not used to constrain the model parameters. 



12 More et al. 




(a) I 



12 13 14 15 

log (M/h-'Mj 



0.15 
0.1 
0.05 - 




1 1 1 1 1 1 1 1 1 



'log M. 



u.v 



(b) 

1 1 1 1 1 



0.1 0.2 0.3 0.4 



00 

o 




12 13 14 15 
log (M/h-'Mj 




9.5 



12 13 14 15 

log (M/h->Mo) 




0.05 



0.1 0.2 0.3 0.4 




12 13 14 15 
log (M/h->Mo) 



Figure 8. Constraints on the model ingredients obtained from the MCMC The 68 and 95 percent confidence intervals are indicated as 
blue and purple regions, respectively. Panel (a) shows the average stellar mass of red centrals as a function of halo mass, while panel (b) 
shows the posterior distribution of the scatter criog jv/^ in this relation. Panels (d) and (e) show the same but for blue centrals. Panel (c) 
shows the constraints on the average number of satellites as a function of halo mass, and finally, panel (f) shows the constraints on the 
fraction of red centrals as a function of halo mass. 



M. < 10^°-^ h-'^MQ, the average halo mass of red centrals 
is virtually identical to that of blue centrals, certainly within 
the errorbars. At larger stellar masses, red centrals reside on 
average in more massive haloes than blue centrals. In con- 
trast, at fixed luminosity (see left-hand panel), the average 
halo mass of red centrals is always systematically larger that 
that of blue centrals, by more than a factor of three. Hence, 
we conclude that stellar mass is a better indicator of halo 
mass than luminosity, especially at the low mass end. How- 
ever, we also stress that there still is a significant amount of 
scatter in the relation between the stellar mass of a central 
galaxy and the mass of the halo in which it resides. 



Finally, recall that we have assumed that the number 
of satellites in a given halo, {Ns)m, is independent of the 
color of its central. In order to test the possible implications 
of this assumption, we have repeated the above analysis al- 
lowing for independent (A's) m for haloes with red and blue 
centrals (each parameterized with Eq. [21]). This adds an 
additional two fi^ee parameters to the model, bringing the 
total to 16. We find that the resulting constraints are per- 
fectly consistent with {Ns)m being identical for red and blue 
centrals, and that all other constraints are similar to what 
we presented above. Hence, we conclude that our assump- 
tion that {Ns)m is independent of the color of its central is 
supported by the data and does not bias our results. 



4.3 Comparison with other studies 

It is interesting and important to compare our constraints 
on the MLR and MSR of central galaxies to those obtained 
using different, independent data sets and methods, includ- 
ing galaxy group catalogues, galaxy-galaxy lensing, galaxy 
clustering, halo abundance matching, and other studies of 
satellite kinematics. In order to enable a fair and mean- 
ingful comparison, whenever required we have adapted the 
results in the literature to match the definitions of halo 
mass and stellar mass used in this paper. In particular, we 
follow ,Hu fc K ravtso v (2003) to convert betwee n different 
definitions of halo mass and use the results of et al.l 
and ^Borch et al.l (|2006l ) to convert stellar masses to 



the 



Kroup a (2001) initial mass function adopted here. 



lYang et al.l ()2007l ) studied galaxy groups in the SDSS, 
which they assigned halo masses based upon either the total 
stellar mass or the total luminosity content of each group. 
We use their group catalogue to investigate the MLR and 
MSR of central galaxies, with and without the split in red 
and blue by colour. The solid and dashed lines in Figs. 1101 
and 1 1 1 1 correspond to the MLR and MSR of central galaxies 
in their group catalog, where the halo masses have been 
assigned using the total stellar mass and total luminosity 
content of the groups, respectively. 

iMandelbaum et al.l (2006) measured the galaxy-galaxy 
lensing signal in the SDSS for galaxies stacked by lumi- 
nosity and by stellar mass. The galaxies were split into 
early and late types based upon their morphology. Here 
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Figure 9. Same as Fig.gl but for the MSR, rather than the MLR. 



we make the crude assumption that this is equivalent to 
our split in red and blue; although clearly an oversimpli- 
fication, we do not believe that this invalidates a compar- 
ison with our results. Modelling their data, they obtained 
the MLR and MSR indicated in Figs. [10] and [TT] by open 
squares with errorbars. The same ex ercise was repeated in 
iMandelbaum. Seliak. fc Hiratal (|2008l '). but by stacking iso- 
lated galaxies. The results of this analysis are shown as filled 
squares with errorbars (95 percent) in Fig. 1111 In the same 
figure, we also show th e galaxy-galaxy lensing results of 
ISchulz. Mandelbaum. fc Padmanabhan (,2009.) as hexagons 
with errorbars (95 percent). 

The abundance and clustering of galaxies holds im- 
portant information regarding the halo occupation statis- 
tics of these galaxies. Using the observed luminosity func- 
tion and clustering properties of galaxies in the SDSS, 
ICacciato et al.l (|2009l ) constrained the conditional luminos- 
ity function, $(L|M), which describes the average number 
of g alaxies of luminosity L that resid e in a halo of mass 
M (|Yang. Mo. fc van den Boschll2003l ). We use their best 
fit parameters to calculate the MLR of central galaxies and 
show it using ope n circlejfl in th e top left-hand panel of 
Fig. 1101 Note that ICacciato et al.l (|2009l ) have also shown 
that their halo occupation model is able to reproduce the 
galaxy-galaxy lensing signal obtained by stacking central 
galaxies based upon their luminosities. 



* ICacciato et al. I ||2009|) also used results from the group cata- 
logue of I Yang et al. to constrain their model, which ex- 

plains why (i) their errorbars are extremely small and (i i ) thei r 
results are in perfect agreement with those of lYang et al.l | |2007| ) 



The relation between galaxy stellar mass and halo mass 
can also be obtained by matching the abundance of galax- 
ies to the abundanc e of haloes and su b haloe s that host 
these galaxies (e.g . Conrov et al. 2006 . 200g| : IGuo et all 
I2OO9I : iMoster et a"lTl2010l : iBehroozi et al.ll2010l V The resuhs 



from this (sub) halo abundance matching technique are often 
quoted as the mean of the distributions P{\ogMt\M). We 
asked the respective authors to provide us log(M)(M«) to 
enable a fair compariso n. In the top left p anel of Fig llll 
we show the res ults of iMoster et al.l l|201Cl|, dotted line), 
IGuo et al] l|2009l . dot-long-dashed line) and IBehroozi et al.l 
l|20ld . dot-dashed hue). 



Fina lly, using data from the SDSS and the DEEP2 sur- 
vey, Co nrov et al.l ^2007) used the kinematics of satellite 
galaxies to determine the evolution of the stellar mass-to 
light ratio of central galaxies from z ^ 1 to 2; ~ 0. They 
measured and modelled the radial dependence of the veloc- 
ity dispersion (in contrast to the aperture averaged velocity 
dispersions used in this paper) to infer the average halo mass 
as a function of the stellar mass of the central galaxy. The 
halo mass— stellar mass relation for all central galaxies thus 
obtained from their analysis at 2 ~ is shown in Fig. [TT] 
using circles with errorbars. 

In all panels of Figs. llOl and fTTI the shaded areas show 
the 95 percent confidence intervals obtained in this paper 
using the kinematics of satellite galaxies. Overall, the re- 
sults obtained using all these very different methods are 
in remarkably good agreement with ea ch other (see also 
IBehroozi et al. I I2OIOI : iDutton et aLll2O10D . and with the re- 
sults obtained here. However, there also are a few discrep- 
ancies, which we discuss below. 




log (Lyh-2 L,,) 



Figure 10. Comparison of our MLR constraints with other constraints from the literature. The shaded regions show the 95 percent 
confid ence intervals that we obtained from the analysis of satellite kinematics. The circles shows the MLR obtained bv lCacciato et aU 
l|2009l l by using galaxy abundance and clustering measurements (the corresponding 95 percent confident intervals are smaller than the 
circles), the squares with errorbars (95 percent confidence intervals) show the MLR obt ained by iM andclba um ct al.. i200& ) using weak 
lensing. The solid and dashed lines show the MLR obtained from the group catalogue of lYang et al.l l|2007l l. 



With regard to the MLR, the galaxy group catalogues 
of lYang et all ||2007|) and the g alaxy-galaxy lensing analy- 
sis of Mandelbaum et all l|2006l ) yield halo masses around 
red centrals that are in good agreement with each other. 
But at L ~ 10^'^h~^ Lq, the halo masses inferred from 
these two methods are roughly a factor two lower compared 
to the results obtained here. The three methods, however, 
agree fairly well at the bright end. Somewhat surprisingly, 
the MLR of all centrals (upper left-hand panel of Fig. llOp 
does not reveal any discrepancy between the masses in- 
ferred from satellite kinematics, versus those inferred from 
either clustering (results of Cacciato ct al. 200^) or galaxy 
group catalogues (results of lYang et all 120071 '). In case of 
the MSR of red centrals, at the low stellar mass end the 
group catalogue results agree with the weak lensing results 
and are again a factor two to three lower than the results 
obtained here. On the other hand, at the bright end, our 
results agree fairly well with the weak lensing results while 
the group catalogue results are roughly larger by a factor 



of 1.6. It is also worth noting that the weak lensing re- 
sults are not fully consistent among each other at the in- 
termediate and low mass ends. In particular, our results 
are in excellent agreement with the weak lensi ng analysis of 
ISchulz. Mandellaaum. fc Padmanabhan (l2009h and th e low 
mass point of Mandelbaum. Seliak. 8z Hiratal (|2008l ). For 
the MSR of all centrals (upper left-hand panel of Fig. Ilip . 
our analysis of satellite kinematics once again yields halo 
masses around low-mass centrals that are ~ 0.3 dex larger 
than those inferred using either subhalo abundance match- 
ing or galaxy group catalo gues. It is notewort hy, though, 
that the results obtained bv lConrov et al.l ()2007l ). which are 
also based on satellite kinematics, are actually in good agree- 
ment with our results. Finally, we note that for blue centrals 
there is no clear indication of any systematic discrepancy, 
except perhaps at the massive end. However, since the mas- 
sive (bright) end of the galaxy mass (luminosity) function is 
dominated by red centrals, the corresponding number statis- 
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Figure 11. Comparison of our MSR constraints with other constraints from the literature literature. The shaded regions show the 95 
percent confidence intervals that we obtained from the an alysis of satellite kinematics. The circles with errorbars (68 percent confidence 
intervals) shows the MSR obtained bv lConrov et al.] ||2007| ) using satellite kine matics. The open squares, t h e solid squares and the hexagons 
with e rrorb ars (95 percent confidence intervals) sh o w the MSR obtained by iMandelbaum et al.l l|2006l ). [Mandelbaum. Seliak. &: Hiratal 
1I2OO8I ') and ISchulz. Mandelbaum, fc PadmanabhanI ||2009^ respectively, using weak lensing. The solid and dashed lines show the MSR 
obtained from the group catalogue of lYan^'eT^ The (mage nta) dotted line, the ( red) dot-dashed line a nd th e (cyan) do t -Iong - 

dashed line show the results from abundance matching studies from [Moster et aL I l l201(f) ,[ Behroozi et al. I | |2010|'| and IGuo et al.l 1I2OO9I ') 
respectively. 



tics are poor resulting in large errorbars. Consequently, we 
do not consider this discrepancy significant. 

To summarize, Figs. [10] and [11] indicate that tremen- 
dous progress has been made in recent years in constraining 
the galaxy-dark matter connection, with different techniques 
yielding MLRs and MSRs that are in fairly good agreement 
with each other, typically within a factor of two. While it is 
difficult to make any robust statement about possible sys- 
tematics, we acknowledge that there is a hint that satellite 
kinematics yields halo masses around low mass centrals that 
are systematically larger than most other methods, espe- 
cially around red centrals. Although we certainly can't rule 
out any systematics in the other methods, we briefly discuss 
a potential problem with satellite kinematics. 

Recently, ISkibba et all (I2OI0I') analyze d the SDSS 
galaxy group catalogue of [Yang et alj (|2007l ) and showed 
that in a significant fraction of groups (ranging from ~ 25 



percent at the low mass end to ~ 40 percent at the massive 
end) the brightest group member is a satellite galaxy rather 
than a central. As discussed at length in their paper, this 
could cause satellite kinematics to overestimate halo masses 
by as much as a factor of ~ 1.6. However, we consider it un- 
likely that this explains the systematic offset between satel- 
lite kinematics and other methods because of the following 
two reasons: First of all, the effect is expected to be largest at 
the massive end, and to be negligible at the low mass end, 
opposite to the trends seen in Figs, [in] and 1111 Secondly, 
the factor ~ 1.6 overestimate only occurs if the probability 
7'bnc(A^) that the brightest (most massive) galaxy in a halo 
of mass M is not the central galaxy is independent of the 
luminosity (stellar mass) of the central galaxy. However, if 
haloes of mass Af in which the central galaxy is fainter (less 
massive) than the average for haloes of that mass are more 
prone to having a brighter (more massive) satellite, which 
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Figure 12. The ratio of the satellite-weighted velocity dispersion, 
(Tsw, measured around hosts that have 7?.i2 (the ratio of their own 
stellar mass to that of their most massive satellite) greater than 
or equal to 7?.JJ,'", to the ratio of the satellite- weighted velocity 
dispersion measured around all hosts {u^^). 



does not se em unreasonable, then the effect can be much 
weaker fsee lSkibba et al.ll2010l . for details). 

To investigate the possibility that our results for red 
centrals at the low stellar mass end are affected by haloes 
in which the central galaxy is not the most massive, we per- 
form the following test. We take all hosts from Sample SR 
that have stellar masses in a bin of width 0.16 dex centered 
around log{M,/ /i~'^Mq) = 10.39. For each host, we deter- 
mine the ratio, 7?.i2, of its stellar mass to that of its most 
massive satellite. Under the assumption that a larger 7?.i2 
implies a larger probability that this host is a true central, 
we proceed as follows. We sort the hosts in decreasing or- 
der of 7?.i2, and create 5 samples by discarding the final 
10, 20, 30, 40 and 50 percent of the hosts (plus their satel- 
lites), respectively. For each of these samples we compute 
the satellite-weighted velocity dispersion of the remaining 
satellites, which we plot in Fig. [12] as a function of the min- 
imum value of TZi2- Clearly, within the errorbars there is 
no indication for a systematic trend between CTsw and 7?. 12". 
Although this does not rule out that our satellite kinematics 
are affected by the possibility that a certain fraction of our 
host galaxies in reality are satellites rather than centrals, it 
certainly makes it less likely. 



5 SUMMARY 

We have used the kinematics of satellite galaxies in order 
to probe the halo mass- luminosity relation (MLR) and the 
halo mass-stellar mass relation (MSR) of central galaxies. 
For this purpose, an iterative selection criterion was first 
applied to select central and satellite galaxies from a volume- 
limited sub-sample of the SDSS. The resulting sample con- 
sist of ~ 6200 satellites around ~ 3900 centrals, making it 
the largest volume-limited sample of central-satellite pairs 
used to date for studies of satellite kinematics. Since the 
number of satellite galaxies around any individual central 
galaxy is small, a stacking procedure was used to combine 
the velocity information of the satellite galaxies of centrals 
with similar luminosities or stellar masses. A detailed mod- 
elling procedure, outlined in Paper I, was then used to infer 
both the average scaling relation between halo mass and the 
central galaxy property and its scatter. 

As expected, and in qualitative agreement with many 
other studies, we find that more luminous (more massive) 
centrals reside in more massive haloes. In addition, we find 



that the MLR of central galaxies is different for central 
galaxies of different colour: red centrals, on average, occupy 
more massive haloes than blue centrals of the same luminos- 
ity. Consequently, the scatter in the MLR of central galaxies 
is at least partly correlated with the colour of the central 
galaxy. 

When stacking central galaxies according to their stel- 
lar masses, we find the difference in the mean MSRs of red 
and blue centrals to be less pronounced than in the case of 
the MLR. In particular, for Af, < lO^" '^ /i"^Mq, the aver- 
age halo masses of red and blue centrals are not significantly 
different. We thus conclude that the stellar mass of a central 
galaxy is a more reliable indicator of its halo mass than its 
(r-band) luminosity. However, even the MSR has a signifi- 
cant amount of scatter of the order of ~ 0.2 dex in stellar 
mass at fixed halo mass. This translates into a scatter in halo 
mass at fixed stellar mass that increases from ~ 0.1 dex at 
M. ~ 4 X 10^ h-^MQ to ~ 0.4 dex at M, ~ 2 x lO" /i'^Mq. 

We compared our constraints on the MLR and MSR 
of central galaxies with a number of other, independent 
constraints coming from the SDSS galaxy group catalogue 
of Yang et al. ( 2007 ) , the g alaxy-galaxy lensing analyses of 
iMandeibaum et all (|2006|), iMandelbaum. Seliak. fc Hiratal 



(|2008l ) and ISchulz. Mandelbauni. fc PadmanabhanI ll2009L 
the galaxy clustering analysis of lCacciato et al.l ll2009r), sub - 
halo abundance matc hing studies of Moster et al. l|2010t ). 
IGuo et al.1 (|2009l ) and iBehroozi et all (|2010|'). and t he anal- 
ysis of satellite kinematics by IConrov et atl l|200'iD . Over- 
all, there is good agreement among all these different stud- 
ies, with a typical study-to-study scatter of less than a 
factor two, which is comparable to the typical 2a errors 
quoted by most of these studies. There is some indica- 
tion, though, that satellite kinematics yield halo masses 
around low-mass centrals (M* <, 3 x 10^" H'^Mq) that are 
systematically higher by a factor 2-3 that most other 
methods, although we emphasize that our results are in 
pe rfect agreement with the galaxy-galaxy lensin g analysis 
of ISchulz. Mandelbaum. fc PadmanabhanI (|2009l ). As dis- 
cussed in the text, we do not believe that our results are 
significantly affected by the fact that not all central galaxies 
are the brightest (or most massive) gal a xies i n their dark 
matter haloes, as shown bv lSkibba et al.l (|2010l ). In fact, de- 
tailed tests with mock galaxy redshift surveys, presented in 
Paper II, have revealed no systematic effects for our method 
of analysis. 

We conclude that the overall level of agreement re- 
garding the MLR and MSR among all different techniques 
indicates that we are converging on an accurate and re- 
liable de scription of the galaxy-da r k matter connection 
(see also van den Bosch et al. I I2OO7I : iBehroozi et al] I2OI0I : 
iDutton et al.ll2010l ). In addition to an overall agreement re- 
garding the means of the MSR and MLR, to well within 
a factor of two, there is also good agreement regarding 
the amount of scatter; as demonstrated in Paper II, the 
scatter of ~ 0.2 dex in luminosity or stellar mass at 
fixed halo mass is in ex cellent agreement with independen t 
constraints obtained by lYang. Mo fc van d en Bosc hl (|2008l ). 
ICacciato et al.l (|2009l ) and Cooray (2006), and with predic- 
ti ons from the s e mi-an alytical model for galaxy formation 
of ICroton et al.l (|2006l ) . This overall level of agreement is 
an admirable achievement, which will prove invaluable not 
only for furthering our understanding of galaxy formation. 
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but also for using galaxies to probe the cosmic density field 
and to constrain cosmological parameters. 



6 ACKNOWLEDGMENTS 

It is a pleasure to thank Eric Bell, Kris Blindert and 
Anupreeta More for useful discussions, and Qi Guo, Peter 
Behroozi and Benjamin Moster for providing results from 
their abundance matching. For a significant duration of this 
project, SM and MC were members of the IMPRS for As- 
tronomy & Cosmic Physics at the University of Heidelberg. 



REFERENCES 

Adelman-McCarthy J.K., et al., 2006, ApJS, 162, 38 
Baldry I.K., Glazebrook K., Brinkmann J., Ivezic Z., Lup- 

ton R.H., Nichol R.C., Szalay A.S., 2004, ApJ, 600, 681 
Becker M.R., et al., 2007, ApJ, 669, 905 
Behroozi, P.S., Conroy, C, & Wechsler, R.H. 2010, preprint 

(arXiv:1001.0015) 
Bell E.F., Mcintosh D.H., Katz N., Weinberg M.D., 2003, 

ApJS, 149, 289 
Blanton M.R., et al., 2003a, ApJ, 592, 819 
Blanton M.R., et al., 2003b, AJ, 125, 2348 
Blanton M.R., Eisenstein D., Hogg D.W., Schlegel D.J., 

Brinkmann J., 2005, ApJ, 629, 143 
Blanton M.R., et al., 2005, AJ, 129, 2562 
Borch A., et al., 2006, A&A, 453, 869 
Brainerd T.G., Specian M.A., 2003, ApJ, 593, L7 
Brown M. J. I., et al., 2008, ApJ, 682, 937 
Cacciato M., van den Bosch F.C., More S., Li R., Mo H.J., 

Yang X., 2009, MNRAS, 394, 929 
Carlberg R.G., Yee H.K.C., EUingson E., Abraham R., 

Gravel P., Morris S., Pritchet C.J., 1996, ApJ, 462, 32 
Carlberg R.G., Yee H.K.C., EUingson E., 1997, ApJ, 478, 

462 

CoUister A.A., Lahav O., 2005, MNRAS, 361, 415 
Conroy C, et al., 2005, ApJ, 635, 982 
Conroy C, Wechsler R.H., Kravtsov A.V., 2006, ApJ, 647, 
201 

Conroy C, et al., 2007, ApJ, 654, 153 
Conroy C, Wechsler R.H., 2009, ApJ, 696, 620 
Cooray A., 2006, MNRAS, 365, 842 
Croton D.J., et al., 2006, MNRAS, 365, 11 
Dai X., Kochanek C. S., Morgan N. D., 2007, ApJ, 658, 
917 

Button A. A., Conroy C, van den Bosch F.C., Prada F., 
More S., 2010, MNRAS, submitted 

Erickson L.K., Gottesman S.T., Hunter J.H., Jr., 1987, Na- 
ture, 325, 779 

Gavazzi R., Treu T., Rhodes J.D., Koopmans L.V.E., 

Bolton A.S., Buries S., Massey R.J., Moustakas L.A., 

2007, ApJ, 667, 176 
Guo, Q., White, S., Li, C, & Boylan-Kolchin, M. 2009, 

preprint (arXiv:d9CI9:4305) 
Hu W., Kravtsov A.V., 2003, ApJ, 584, 702 
Kravtsov, A.V., Berlind, A. A., Wechsler, R.H., Klypin, 

A.A., Gottlober, S., AUgood, B., & Primack, J.R. 2004, 

ApJ, 609, 35 
Kroupa P., 2001, MNRAS, 322, 231 



Li C, Kauffmann G., Jing Y.P., White S.D.M., Borner G., 

Cheng F.Z., 2006, MNRAS, 368, 21 
Li R., Mo H. J., Fan Z., Cacciato M., van den Bosch F. C, 

Yang X., More S., 2009, MNRAS, 394, 1016 
Maccio A. v., Button A. A., van den Bosch F.C., Moore B., 

Potter D., Stadel J., 2007, MNRA S, 378, 55 

McKay T.A., et al., 2001, preprint, ||astr o-ph/0108013| | 

McKay T.A., et al., 2002, ApJ, 571, L85 

Mandelbaum R., Seljak U., Kauffmann G., Hirata CM., 

Brinkmann J., 2006, MNRAS, 368, 715 
Mandelbaum R., Seljak U., Hirata C. M., 2008, JCAP, 8, 

6 

More S., van den Bosch F.C., Cacciato M., Mo H.J., Yang 

X., Li R., 2009, MNRAS, 392, 801 
More S., van den Bosch E.G., Cacciato M., 2009, MNRAS, 

392, 917 

More S., 2009, Ph. D. thesis, Ruprechts-Karls-Universitat, 
Heidelberg 

Moster B.P., Somerville R.S., Maulbetsch C, van den 
Bosch F.C., Maccio' A.V., Naab T., Oser L., 2010, ApJ, 
710, 903 

Navarro J.F., Frenk C.S., White S.D.M., 1997, ApJ, 490, 
493 

Norberg P., Frenk C.S., Cole S., 2008, MNRAS, 383, 646 
Parker L.C., Hoekstra H., Hudson M.J., van Waerbeke L., 

Mellier Y., 2007, ApJ, 669, 21 
Prada F., et al., 2003, ApJ, 598, 260 

Rubin V.C., Ford W.K., Jr., Thonnard N., Burstein D., 

1982, ApJ, 261, 439 
Rykoff E.S., et al., 2008, ApJ, 675, 1106 
Schulz A. E., Mandelbaum R., Padmanabhan N., 2009, 

preprint (arXiv:0911.2260 ) 
Seljak U., 2000, MNRAS, 318, 203 

Skibba R., Sheth R. K., Connolly A. J., Scranton R., 2006, 

MNRAS, 369, 68 
Skibba R. A., Sheth R. K., 2009, MNRAS, 392, 1080 
Skibba, R.A., van den Bosch, F.C., Yang, X., More, S., 

Mo, H., & Fontanot, F. 2010, MNRAS submitted, preprint 

(arXiv: 1001. 4533) 
Spergel D.N., et al., 2007, ApJS, 170, 377 
Tinker J.L., Weinberg D.H., Zheng Z., Zehavi I., 2005, ApJ, 

631, 41 

Tinker J., Kravtsov A.V., Klypin A., Abazajian K., Warren 
M., Yepes G., Gottlober S., Holz D.E., 2008, ApJ, 688, 709 

van den Bosch F.C., Yang X., Mo H.J., 2003, MNRAS, 340, 
771 

van den Bosch E.G., Norberg P., Mo H.J., Yang X., 2004, 

MNRAS, 352, 1302 
van den Bosch F.C., et al., 2007, MNRAS, 376, 841 
van den Bosch F.C., Aquino D., Yang X., Mo H.J., Pasquali 
A., Mcintosh D.H., Weinmann S.M., Kang X., 2008, MN- 
RAS, 387, 79 

Yang X., Mo H.J., van den Bosch F.C., 2003, MNRAS, 339, 
1057 

Yang X., Mo H.J., van den Bosch E.G., Weinmann S.M., 

Li C, Jing Y.P., 2005, MNRAS, 362, 711 
Yang X., Mo H.J., van den Bosch E.G., Pasquali A., Li C, 

Barden M., 2007, ApJ , 671, 153 
Yang X., Mo H.J., van den Bosch E.G., 2008, ApJ, 676, 248 
Yang X., Mo H.J., van den Bosch E.G., 2009, ApJ, 695, 900 
Yoo J., Tinker J.L., Weinberg D.H., Zheng Z., Katz N., 

Dave R., 2006, ApJ, 652, 26 



18 More et al. 



York D.G., et al., 2000, AJ, 120, 1579 
Zaritsky D., Smith R., Frenk C, White S.D.M., 1993, ApJ, 
405, 464 

Zaritsky D., White S.D.M., 1994, ApJ, 435, 599 
Zaritsky D., Smith R., Frenk C, White S.D.M., 1997, ApJ, 
478, 39 

Zehavi I., et al., 2004, ApJ, 608, 16 
Zehavi I., et al., 2005, ApJ, 630, 1 



APPENDIX A: SAMPLE SELECTION 
FUNCTION 

Here we illustrate how our use of a volume limited sam- 
ple of galaxies, complete in luminosity, leads to a sample 
that is incomplete in stellar mass. Our aim is to character- 
ize the sample selection function, S{M,, C), which describes 
the fraction of galaxies of stellar mass M* and colour C that 
make it into our sample. This sample selection function is 
required for modelling the satellite kinematics as function of 
the stellar mass of their host galaxy (see Section [3. 3. 2|) 

The left-hand panel of Fig. lAll shows the distribution 
of galaxies in the "'^r-band luminosity-redshift plane. Red 
and blue galaxies are indicated by red and blue dots, respec- 
tively. The apparent magnitude limit, mj!™ — 17.77, of the 
spectroscopic sample results in an absolute magnitude limit 
given by 

"■^M"'" - 51og/i = 17.77- DM(z) - ko.i{z) 

+1.62(2-0.1). (Al) 

Here ko,i(z) is the k-correction to r edshift z = 0.1 , 1.62 is 
the evolution correction factor from iBlanton et al.l (|2003al ). 
and DM(2:) is given by 

DM(z) = 51ogDL(z) + 25.0, (A2) 

where Dl{z) is the luminosity distance of the galaxy in 
/i~^Mpc. The redshift depende nce of the k-correcti ons is 
fairly well reproduced by ( van den Bosch et al. II2OO8I ) 

fco.i(2) = 2.51og(i±^) . (A3) 

The above equations yield the solid black line in the left- 
hand panel of Fig. lAll Note that a very small fraction of 
galaxies fall below this predicted limit. This is because the 
k-correction also depends on the colour of the central galaxy 
which we have not accounted for. However this effect can 
safely be ignored for the purpose of this paper. 

Galaxies with z ^ 0.072 and °-^Lr > lO'''^ (i.e., 
°-^M"™ - 5 log /i ^ -18.99) make up the volume hmited 
sample used in this paper (indicated by dashed lines in 
the left-hand panel of Fig. IA1[) . The right-hand panel of 
Fig. lAll shows the distribution of galaxies in this volume lim- 
ited sample in the stellar mass-redshift plane. Note that the 
sharp cut in luminosity translates into a colour-dependent 
cut in stellar mass; bluer galaxies have a lower cut-off in 
stellar mass, causing the low-mass end of the sample to be 
completely dominated by blue galaxies. 

We remind the reader that in our analysis, galaxies were 
assigned stellar masses using the "'^{g — r) colour and the 
r-band magnitude using Eq. (dl). The ° V band magnitude 
limit we have used for our volume limited sample translates 
into a "'"r-band limit given by 
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Figure A2. The colour-stellar mass diagram of galaxies in our 
volume-limited sample. The black solid line shows the boundary 
that we use to split the galaxy population in 'red' and 'blue' 
galaxies. The green dashed line shows the analytical prediction 
(Eq. IA6n for the selection effect that results due to our use of 
a volume-limited sample complete in luminosity. The resulting 
selection functions, 5(C, A/»), for red and blue galaxies are in- 
dicated by the (red) dashed and (blue) dot-dashed lines, respec- 
tively. 

O.O^lim ^0.1 ^^lim ^ j^^^^^^^ _ ^^^^^j _ Q ,^g2 , (A4) 

where the term in square brackets is the difference in k- 
corrections between redshift 0.0 and 0.1, and the last term 
is the difference in the evolution corrections between these 
redshifts. The k-co rrections to redshift 0.0 can be reasonably 
well described by (|van den Bosch et al. II2OO8I ) 

ko.(){z) = 2.51og(l + z) + 1.5 z [°-°{g - r) - 0.66] (A5) 

Combining Eqs. (|X4)) . (|X5)) . (|X3|l with Eq. ([TJ gives that, at 
fixed " "(g — r) colour, the stellar mass limit in the volume 
limited sample varies with " "(gr — r) and redshift as 

\og{M'r) = 8.9812 - log [y4^t^ 
[ 1.1 (1 + z) 

-0.396 z + (1.097 + 0.6 z f-°{g - r) (A6) 

In Fig |A2l we plot the distribution of galaxies in the 
colour-stellar mass diagram. The colour cut which was used 
to assign colours to our sample of galaxies is shown by the 
black solid line (see Eq. [T]). The green dashed line show the 
limit expressed in Eq. (|A6|) assuming z = 0.072. At fixed 
stellar mass, only a fraction of galaxies (those that lie below 
the green line), are part of the volume limited sample. More 
than half of the red galaxies drop out of the sample at stellar 
masses below lO^'® /i'^Mq. 

To calculate the selection function S{M,,C), we use 
the entire flux limited catalogue of galaxies to populate the 
"■"(g — r)— stellar mass plane, this time assigning each galaxy 
a weight equal to one over the maximum volume to which 
this galax;y could be seen given the r-band flux limit of 17.77. 
Using small bins in stellar mass, we calculate S{M,,C) in 
each bin by dividing the sum of weights of galaxies of a 
particular colour that lie below log(Af''™) with the sum of 
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Figure Al. The left panel shows the distribution of SDSS galaxies in the luminosity-redshift plane. Red and blue galaxies are indicated 
by red and blue dots, respectively. The volume-limited sample of galaxies used for our analysis of satellite kinematics is enclosed by the 
dashed lines. The right hand panel shows the distribution of galaxies in this volume-limited sample in the stellar mass-redshift plane. 



weights of all galaxies of that colour in the bin under con- 
sideration. The resulting selection functions for red and blue 
galaxies are plotted in Fig. IA2l as red, long-dashed and blue 
dot-dashed lines, respectively. 



